Theory of incompressible MHD turbulence with scale-dependent alignment and cross-helicity 
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A phenomenological anisotropic theory of MHD turbulence with nonvanishing cross-helicity is constructed 
based on Boldyrev's phenomenology and probabilities p and q for fluctuations 5v±_ and 6b± to be positively or 
negatively aligned. The positively aligned fluctuations occupy a fractional volume p and the negatively aligned 
fluctuations occupy a fractional volume q. Guided by observations suggesting that the normalized cross-helicity 
o c and the probabilities p and q are approximately scale invariant in the inertial range, a generalization of 
Boldyrev's theory is derived that depends on the three ratios w + /w~ , e + /e~, and p/q. It is assumed that 
the cascade process for positively and negatively aligned fluctuations are both in a state of critical balance and 
that the eddy geometries are scale invariant. The theory reduces to Boldyrev's original theory when a c = 0, 
e + = e~, and p — q and extends the theory of Perez and Boldyrev when a c 7^ 0. The theory is also an 
anisotropic generalization of the theory of Dobrowolny, Mangeney, and Veltri. 

PACS numbers: 



I. INTRODUCTION 

Phenomenological theories of incompressible MHD turbu- 
lence that take into account the anisotropy of the fluctuations 
with respect to the direction of the mean magnetic field Bo 
were pioneered by Goldreich & Sridhar and others in the 
1990s. The influencial and somewhat controversial theory 
of Goldreich & Sridhar 01 0] established the idea that the 
timescale or coherence time for motions of a turbulent eddy 
parallel and perpendicular to Bo must be equal to each other 
and this unique timescale then defines the energy cascade 
time. This concept, called critical balance, leads immediately 
to the perpendicular energy spectrum E(k±) oc k ± 5 ^ 3 and to 



the scaling relation fen esc & 
turbulent eddies. 
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describing the anisotropy of the 



The decade following the publication of the paper by Gol- 
dreich & Sridhar in 1995 was a time when significant ad- 
vances in computing power were brought to bear on com- 
putational studies of MHD turbulence. Simulations of in- 
compressible MHD turbulence during this time showed that 
when the mean magnetic field is strong, Bq {8b) 2 , the 
perpendicular energy spectrum exhibits a power-law scaling 



closer to k^ 2 than fc ± J/J |3|-|5|]. Motivated by this result, 
Boldyrev modified the Goldreich & Sridhar theory to explain 
the k^ 2 power-law seen in simulations 10,0]. A new con- 
cept that emerged in Boldyrev's theory is the scale dependent 
alignment of velocity and magnetic field fluctuations whereby 
the angle 9 formed by 8v± and 5b± scales like X 1 / 4 in the 
inertial range. This alignment effect weakens the nonlin- 
ear interactions and yields the perpendicular energy spectrum 
E(kj_) oc k ± 3 ^ 2 . Evidence for this alignment effect has been 
found in numerical simulations of incompressible MHD tur- 
bulence it Ht] and in studies of solar wind data ITolfTlll . 
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The phenomenological theory of Galtier et al. 11211 can also 
be used to explain the observed k ± energy spectrum. Us- 
ing a slightly modified critical balance relation that retains 

2/3 

I oc fej_' scaling of the Goldreich & Sridhar theory, 



the k 
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their model admits the k ± energy spectrum (as well as 
other solutions). However, the theory of Galtier et al. lfl2tl 
does not include the scale dependent alignment that arises in 
Boldyrev's theory and, more importantly, is seen in the solar 
wind OSlH. 

The theories discussed so far 01 0, OH] all assume 
that the cross-helicity vanishes and, therefore, these theories 
cannot be applied to solar wind turbulence. When the cross- 
helicity of the turbulence is nonzero it is necessary to take 
into account the cascades of both energy and cross-helicity. A 
generalization of the Goldreich & Sridhar theory to turbulence 
with nonvanishing cross-helicity, also called imbalanced tur- 
bulence, has been developed by Lithwick, Goldreich, & Srid- 
har OH- Other theories of imbalanced turbulence have been 
derived by Beresnyak & Lazarian lfl4ll and Chandran OHO^l- 
However, none of these theories contains the scale dependent 
alignment of velocity and magnetic field fluctuations seen in 
the solar wind. Therefore, to develop a theory that may be 
applicable to solar wind turbulence it is of interest to general- 
ize Boldyrev's theory to incompressible MHD turbulence with 
non- vanishing cross-helicity. An extension of Boldyrev's the- 
ory to imbalanced turbulence has been discussed by Perez & 
Boldyrev 03] • 

The purpose of the present paper is to develop 
a theoretical framework which generalizes the results of Perez 
& Boldyrev and is consistent with solar wind observations. 

Our theory is founded, in part, on two new solar wind ob- 
servations presented in this paper. The first is the observa- 
tion that the normalized cross-helicity cr c , the ratio of cross- 
helicity to energy, is scale invariant in the inertial range. The 
second is the observation that the probabilities p and q that 
fluctuations are positively or negatively aligned, respectively, 
are also scale invariant, that is, these quantities are approxi- 
mately constant in the inertial range. Experimental evidence 
for the scale invariance of a c comes from solar wind observa- 
tions by Marsch and Tu 11811 an d also Figure 1 below, and from 
numerical simulations ltl7l[l9ll . Evidence for the scale invari- 
ance of p and q is shown in Figure 2 below. Assuming these 
quantities are all scale invariant we deduce expressions for the 
energy cascade rates and the rms fluctuations that generalize 
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the results in |01 and |fl7ll and are consistent with the concept 
of scale dependent alignment of velocity and magnetic field 
fluctuations, a concept neglected in other phenomenological 
theories B12l4l5ll . The resulting theory, which is founded on 
the concept of scale-invariance and grounded in solar wind 
observations, contains the theories of Boldyrev J3l and Perez 
and Boldyrev 111711 as special cases, but opens up a broader 
range of physical possibilities. 

Consistent with numerical simulations and solar wind ob- 
servations, in our approach the fluctuations at a given point 
may assume one of two possible states referred to as positively 
aligned Sv± ■ Sb± > and negatively aligned Svj_ ■ 5b± < 0. 
Each state is characterized by its own rms energy v 2 , align- 
ment angle 0, and nonlinear timescale r. Positively aligned 
fluctuations have a characteristic spatial gradient which deter- 
mines their nonlinear timescale and negatively aligned fluctu- 
ations have a different spatial gradient which determines their 
nonlinear timescale. These timescales are estimated from the 
nonlinear terms in the MHD equations as described in sections 
2 and 3. 

Section 2 describes the geometries of velocity and magnetic 
field fluctuations that are either aligned 'f ' or anti-aligned 'J,' 
and these are used to form estimates of the nonlinear terms 
in the MHD equations. From this foundation, estimates of 
the energy cascade times are constructed in section 3 and the 
theory of the energy cascade process is developed in section 
4. The summary and conclusions are presented in section 5. 



II. FLUCTUATIONS IN IMBALANCED TURBULENCE 

Consider velocity and magnetic field fluctuations measured 
between two points separated by a distance A^ in the field 
perpendicular plane. Let v and b denote the fluctuations in the 
plane perpendicular to the local mean magnetic field, where v 
and b are both measured in velocity units. Suppose that v and 
b are aligned with some small angle > and assume, as for 
Alfven waves, that \v\ = \b\. Then w + = v + b is nearly 
aligned with v and w = v — b is nearly perpendicular to v 
as sketched in Figure 1 . It follows from the identity w ± xv ~ 



=pv x b that 



(1) 



where + is the angle formed by w + and v, 8 is the angle 
formed by w~ and v, and is the angle formed by v and 
b. In addition, 0+ + 0~ = tt/2. Following Boldyrev 0], 
suppose that the gradient of the fluctuations is in the direction 
perpendicular to v. In this case, 



and 



(w~ ■ V) 



0+ • v) 



w sm a v sin f 



) + sin 8 + v sin 



Ax 



Ai 



(2) 



(3) 



The time rate of change caused by nonlinear interactions is 
estimated from the relations 



d \w 



+ |2 



and 



dt 



d \w 



dt 



w ■ (w + ■ V)w 



(4) 



(5) 



If v and b are aligned with some small angle 8, then the 
fluctuations are called "positively aligned" and denoted by 'f 
(Figure la). Similarly, if v and —6 are aligned with some 
small angle 8, then the fluctuations are called "negatively 
aligned" or "anti-aligned" and denoted by 'J,' (Figure lb). For 
positively aligned fluctuations, equations d2J — (f5j imply 



dt 2 ~ XT = XT 



and 



d (w-j- ) 2 (w^ ) 2 w^ (w t ) 2 v-f sin0-|- 



dt 



A 



A j 



(6) 



(7) 



where 6f is the angle formed by v and b and quantities with 
the subscript j" describe positively aligned fluctuations. It is 
clear from the middle term in equation © that the time rate of 
change of depends on wZ, consistent with the nonlinear 
terms in the MHD equations, although this dependence is not 
immediately apparent in the last term in ©. For negatively 
aligned fluctuations, equations (O-© imply 



d (wp 2 (wp 2 w l sin 6^ (wj) 2 v±sm8i 



dt 



and 



d (w i ) 2 (it^ ) 2 w^ sin 6^ _ {w i ) 2 vi sin0; 



dt 



A 



A 



(8) 



(9) 



where 9± is the angle formed by v and —b and quantities with 
the subscript J, describe negatively aligned fluctuations. Here, 
< t < tt/2 and < 8 i < tt/2. 

In general, the fluctuations v and b observed at any point 
(x, t) are either positively or negatively aligned. For a point 
(x,i) picked at random, let p and q be the probabilities the 
alignment is positive or negative, respectively {p + q = 1). 
Then, on average, 

d (w + ) 2 1 

~ — [p(w+)\ sin 8 t + q(w+) 2 Vi sin 8^ (10) 



and 

~ -^[p(wf ) 2 w t sin0 t + q(w7) 2 Vl 

where the rms values are defined by 

(w+) 2 =p(w+) 2 +q(w+) 2 , 

-.-\2 „/...--\2 , „/...- \2 



(w y = P (w t y + q(w±y 



sm 8i\ , 
(11) 

(12) 
(13) 
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FIG. 1 : Geometry of the fluctuation vectors v and b for positively aligned fluctuations (a) denoted by 'f and negatively aligned fluctuations 
(b) denoted by '4-'. The gradient is perpendicular to the velocity fluctuation v. The magnitude of v for positively and negatively aligned 
fluctuations are v-f and t>i, respectively. The angles formed by v and b, 8^ and 9i, are both assumed to be small. 



The following relations also hold. For a positively aligned 
fluctuation, assuming = |6|, 



wf ■ wf = 2vi(l + cos0-|-), 



■j^ • UU^ — 6Uf\ 

2/ 

- u/^. — *Vf\ 



(14) 
(15) 



= 2v 2 (l — cos0f), 

and wfw^ = 2u| sin#-|-. The energy of a positively aligned 
fluctuation is v 2 . For a negatively aligned fluctuation 



w 
w 



f-wf = 2vf(l- cos 6{), 
■~ -w^ = 2vf(l + cos 04.), 



(16) 
(17) 



and w^w^ = 2i>| sin The energy of a negatively aligned 
fluctuation is vf. Thus, the rms values (fT2b and (fT3l are 

(w+) 2 = 2[pv%(l + cos6 t )+qvl{l-cos6 l )], (18) 
(w~) 2 = 2[pw|(l -cos0 t ) + cos 04.)]. (19) 

If the angles are small, ^ C 1 and 0j, <C 1, then the small 
parameter 8 can be used to order the terms in equations ( fT8l 
and ( fT9b so that to leading order 

(w+) 2 ~4u|p and (tIT) 2 ~ 4u|g, (20) 

where = 1. This may be derived as follows. Inequations 
( fT8l and ( fl9l ) assume that the angles are both small and then 
substitute 1 + cos 9 ~ 2 and 1 — cos 6* = 2 sin 2 (9/2) to obtain 



and 



(™+) 2 ~4[pv 2 + qv 2 S m 2 (9j2)] (21) 



(w") 2 ~ 4 [pv\ sin 2 (0 t /2) + ^ 2 ] . (22) 



As A_l — > 0, both 6f ->• and #4. — » and, therefore, to 
first order, the terms proportional to sin 2 (9) may be neglected. 



Alternatively, note that 

1 (pv 2 /qv 2 ) + S in 2 (9 l /2) 

w 



(^ 2 M 2 )sin 2 (0 t /2) + l 



(23) 



As discussed below, solar wind observations show that this 
quantity is approximately constant in the inertial range. Now, 
as Of — > and 9^ — > the only way that this can remain 
constant is if pv 2 /qv 2 is bounded away from zero and 



w 



(24) 



This justifies the approximation in Eqn I 

Equation ( f20b shows that at a given scale A_l the total en- 



ergy [( 



,7,+ V 



+ (w ) ]/4 is partitioned into two parts, the en- 



ergy (w + ) 2 /A associated with positive alignment and the en 
ergy (w)~) 2 /4 associated with negative alignment. The nor- 
malized cross-helicity a c is defined as the ratio of the cross- 
helicity to the energy at a given scale and can be written 

_ (^+) 2 - (^-) 2 .... 

For small angles, equations (TTOb and (fTTl i become, to lead- 
ing order, 

d [w+ f 4p^0 t 



dt 2 
d (w~) 2 
dt 2 



4qv 3 



A 1 



(26) 



(27) 



To express these in terms of the rms values w ± , eliminate 
and vi using equation d20l i. This yields 

d (w+) 2 ^ 



dt 2 

8 (w-) 2 
dt 2 



2A ±P V2 ' 
2X ± q 1 / 2 ' 



(29) 
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These estimates shall be used to derive the cascade times. 



III. ENERGY CASCADE TIME 

When nonlinear interactions are strong and a large number 
of Fourier modes are excited, fluctuations occur continuously 
in time and space. During a time r the fractional change in 
the quantity (u> + ) 2 is, from i 



X + (r) 



2r 



(w+) 3 6> T 



W + 9fT 



T <T^ 



(30) 



where r + is the cascade time at the lengthscale Aj_ and the 
tildes have been dropped. Similarly, the fractional change in 
the quantity (il>~) 2 is, from i 



X-(r) 



{vrfe i 

2X±q 



•It 



1/2 



(w 



A,V/ 2 ' 



T < T 



(31) 



where r is the cascade time of w and the tildes have been 
dropped for brevity. Hereafter, the tildes will be omitted and 



+ and w will always represent the rms values. 



According to the definition of the energy cascade time, the 
fractional change \ + is of order unity when the interaction 
time t is equal to the cascade time t + . Therefore, the relations 
([Hi and (ED imply 



A iP 1/2 



W~ l ~9-[ w t>i 

By similar reasoning, equations (O and (O imply 



A 



A 



n 



(32) 



(33) 



Moreover, equations ( f32b . ( 1331 ), and (f20b imply r + = t-j- and 
t~ = tj.. Thus, the energy cascade times for the rms Elsasser 
amplitudes are equal to the energy cascade times for the posi- 
tively and negatively aligned fluctuations. 

For balanced turbulence, <j c — > 0, w + /w~ — s- 1, p = q, 
Of = 9±, and the energy cascade times (1321 reduce to the cas- 
cade time in Boldyrev's original theory |7[]. For imbalanced 
turbulence, a c ^ 0, the cascade times (l32l are different from 
the cascade times ~ X±/w T 9 T in the theory of Perez & 
Boldyrev lTl7ll . The theory presented here is different from the 
theory of Perez & Boldyrev 11711 because the latter theory does 
not take into account the existence of two separate types of 
fluctuations, positively and negatively aligned, with separate 
probabilities of occurrence p and q. Taking this into account 
and also the definitions of the rms amplitudes dT2b and ( TOl ). 
it follows from the preceding analysis that the timescales for 
the rms amplitudes take the form ( 1321 ). 

As pointed out by Kraichnan [20], Dobrowolny, Mangeney, 
and Veltri 112 ill , and others, the energy cascade in MHD turbu- 
lence occurs through collisions between Alfven wavepackets 
propagating in opposite directions along the mean magnetic 
field. In other words, it is the interaction between w + and 
w~ waves that causes the energy to cascade to smaller scales 



in MHD turbulence. Consequently, the cascade time for w + , 
say, should depend on w~ . While it may appear from equa- 
tions (f32t— <f33T> that the timescale for w + fluctuations depends 
only on w + and, therefore, the interaction with the w~ waves 
is absent, this is not true. The interactions are still present 
in the expressions ( [321 and d33l through the dependence on 
the angles and other parameters as will be shown in the next 
section. 



IV. THEORY OF THE ENERGY CASCADE PROCESS 

Assuming there is no direct injection of energy or cross- 
helicity within the inertial range and there is no dissipation 
of energy or cross-helicity within the inertial range, the en- 
ergy cascade rate e and the cross-helicity cascade rate e c are 
scale-invariant in the inertial range. It follows that the energy 
cascade rates for the two Elsasser variables ± e c are 

also scale-invariant. The theory of the energy cascade process 
is based on Kolmogorov's relations 



(w+) 2 
~2t+~ 



and 



2t- 



(34) 



where the non-zero constants e + and e~ are the energy cas- 
cade rates per unit mass for the two Elsasser variables w + and 
w~, respectively. These equations describe the conservation 
of energy flux in fc-space (Fourier space). In addition to Kol- 
mogorov's relations, there are two observational constraints 
that must be taken into account. 

Solar wind observations show that the energy and cross- 
helicity spectra of the turbulence follow approximately the 
same power law in the inertial range (Fig. 2) which implies 
that the normalized cross-helicity a c is approximately con- 
stant. In other words, the quantity a c is approximately scale 
invariant. Similar results have been found in simulations of 
incompressible MHD turbulence fTil [T7L [l9tl . In particular, 
the 3D simulations of Perez & Boldyrev [ 17] indicate that the 
perpendicular Elsasser spectra are proportional to each other 
in Fourier space. Solar wind observations also suggest that 
the probabilities p and q are approximately scale invariant as 
shown in Fig. 3. These observations will now be taken into 
account in the theory. 

Assuming er c and p are both scale invariant quantities, then 
w + /w~, Wf/uj., t + /t _ , and 9^/6^ are scale invariant by 
equations ( |25] l, d24l . d34l , and d32l , respectively. In all, there 
are six different scale invariant ratios in the theory 



p Vf 

q vi 



(35) 



At most, only three of these are independent, say, the first 
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FIG. 2: Typical energy E and cross-helicity spectra C (trace spec- 
tra) obtained using 3-second plasma velocity and magnetic field data 
from the Wind spacecraft near the orbit of the earth at 1 AU. (a) An 
interval of highly Alfvenic high-speed wind from 3 Jan 1995 09:00 
to 8 Jan 1995 00:00, 4.625 days, (b) A weak high-speed stream em- 
bedded in low-speed wind; 24 Jul 1996 12:00 to 7 Aug 1996 00:00, 
14 days, (c) The normalized cross-helicity a c — C/E as a function 
of frequency. The rapid change in a c near the Nyquist frequency is 
at least partly caused by the FFT processing techniques and may not 
be a real physical effect. 
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FIG. 3: The probabilities p and q obtained from solar wind data by 
integrating the observed probability density function for the angle 
9 from to 7r/2 and from n/2 to ir, respectively. The data was 
acquired by the Wind spacecraft between 8 Jan 1997 and 9 June 1997 
and analyzed using the techniques described in jlljl. Examples of the 
probability density functions can be found in 1 1 ill . 



three. Equations ( 124b . ( 1341 ), and ( f32b imply 




(36) 
(37) 
(38) 



Therefore, the six scale invariant ratios d35l l can all be ex- 
pressed in terms of the first three w + /w~, e + je~ , and p/q. 



To be able to solve Kolmogorov's relations d34i > for w ± it is 
necessary to express the alignment angle Of in terms of w . In 
general, Of can depend on w + , w~, the Alfven speed va, the 
lengthscale Aj_, the cascade rates e + and and the probabil- 
ities p and q. By dimensional analysis, Of must be a function 
of the following six dimensionless quantities 



VA 



w 

VA 



A, 



P, 



(39) 



Moreover, Of must change into Of when w + , e + , and p are in- 
terchanged with w~, e~, and q, respectively, to be consistent 
with the nonlinear terms d28l ) and d29l ). For a theory com- 
posed of power law functions, the only forms that satisfy all 
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these requirements are 

9-f oc 



+ \ a / —\P/ +\ \ 7 



e+A, We- A 



.4 



,1 



Kg", (40) 



qfjf, (41) 



where a, /3, 7, 5, /i, and ^ are constants that must be deter- 
mined by the theory. In addition, there is a leading coefficient 
which is omitted. 

The substitution of (l40l > and ( 1411 into equation d38l yields 
= a + 3, £ = 7 — 1, and v = — 1/2. The parameters are 
further constrained by considering the geometry of the "tur- 
bulent eddies" associated with the fluctuations Vf and uj.. The 
parallel correlation length is defined by Aj = Vat^ and the 
correlation length in the direction of the velocity fluctuation 
is ^ = VfTf. Similarly, the correlation lengths for negatively 
aligned fluctuations are A| = i^tj, and £j_ = v^t^. In the 
plane perpendicular to the local mean magnetic field £ is par- 
allel to v, the gradient direction is perpendicular to v with 
lengthscale Aj_, and the eddy dimensions are £ x Aj_. The di- 
mension parallel to the mean magnetic field is A11. Hence, in 
physical space the turbulent eddies can be visualized as three- 
dimensional structures with dimensions A_i_ X £ X A11. 

The coherence times for longitudinal and transverse mo- 
tions of the eddy must be equal to each other and also to the 
cascade time. This is the critical balance condition of Gol- 
dreich and Sridhar which is also implicit in the work of Hig- 
don ll22tl . Equation d33l and the definitions of the correlation 
lengths in the last paragraph immediately yield the critical bal- 
ance condition 



A 



t 



it 



A 



(42) 



with a similar condition for the negatively aligned fluctuations 



A 



n 



A 1 



va vi 2v±9i 



(43) 



Now consider the eddy geometry. When the mean magnetic 
field is strong enough that /va < 1, then A|| > £ > A^ 
and the eddies are elongated in the parallel direction. The 
condition /va < 1 is assumed hereafter. Equation (l42l 
shows that the aspect ratio in the field perpendicular plane is 
<fo = Aj_/£t = anc ' trie as P ect ra ti° m the parallel direc- 
tion is, from equations (l42l and (|20l ). 



■0t 



k 

a! 



VA 



w ' 



2v A p 1/2 ' 



(44) 



The two aspect ratios will scale in the same way if cf>^ / i/j-f is 
scale invariant. This implies that a = — 1 and 7 = 1/2. The 
assumption that the ratio 4>f/ipf is scale invariant is differ- 
ent from Boldyrev's original approach in which he assumed 
that the alignment angles in and out of the field perpendic- 
ular plane are simultaneously minimized. Nevertheless, our 



assumption retains the spirit of Boldyrev's original theory 
which implies the geometry of turbulent fluctuations are scale- 
invariant. 

Solving Kolmogorov's relation ( f34b using (l32b , d40i >. d4ll . 
and the parameter values obtained so far, one finds 



w 

VA 



±1/2 / £ -x ' 



E ± Xj 



1/4 



iPl) 



-v/A 



and the total energy cascade rate e = (e + + e )/2 is 

(w + w~) 2 
AvaX±_ 

The total energy at scale A^ is 

(w + ) 2 + {w~) 2 W + W~ f W + | n ^ _ 2 




(45) 



(46) 



(47) 



Therefore, the energy cascade rate can be written 




(pq) v - (48) 



Assuming the rms energy v 2 at scale A^ is held constant, the 
terms on the right-hand side describe the dependence of the 
energy cascade rate on the ratios e + /e~ and w + /w~. The 
value of v may be determined by comparison with experiment 
or possibly by further physical considerations. This parame- 
ter does not affect the inertial range scaling laws and is left 
undetermined for the moment. 

At this point it is of interest to return to the expressions 
(|32| > for the cascade times and ask: How do the cascade times 
depend on the rms Elsasser amplitudes? Using the parameter 
values obtained previously, equation (l40b becomes 



(49) 



and the substitution of this result into equation d32t yields 



(50) 



A similar expression holds for t~ so that the ratio t + /t~ 
satisfies ( fTTb . Ignoring scale invariant factors, the preceding 
equation shows that 



(w-) 2 



and 



X±va 
(w+) 2 



(51) 



In this form, the angle dependence has been eliminated. Note 
that the simple estimate r + ~ \±/w~ suggested by the non- 
linear term in the MHD equations is modified by the factor 
va/w~ which accounts for the weakening of nonlinear inter- 
actions caused by scale dependent alignment. The presence of 
this algebraic factor is one of the hallmarks of Boldyrev's orig- 
inal (2006) theory which is generalized here to imbalanced 



7 



turbulence. Remarkably, the relations (l5TT i are identical to 
those in the isotropic theory of imbalanced turbulence devel- 
oped by Dobrowolny, Mangeney, and Veltri; see equation (10) 
in l2lll . Recall that Dobrowolny, Mangeney, and Veltri con- 
cluded from their expressions for the cascade times that steady 
state turbulence with nonvanishing cross-helicity is impossi- 
ble. On the contrary, the theory presented here allows such a 
steady state because the additional coefficients shown in < f50b 
but not dSTb maintain the relation (f37b even when e + ^ e~ . 
Thus, the theory presented here is also a generalization of the 
theory of Dobrowolny, Mangeney, and Veltri 112 ill . 

A remark about the timescales in the theory should be men- 
tioned. If w + > w~, then equation d37l i implies it is possible 
that r + < r~ since there is nothing in the theory that pre- 
vents this. That is, the energy of the more energetic Elsasser 
species may be transferred to smaller scales in less time than 
the energy of the less energetic Elsasser species. This is not in- 
consistent with dynamic alignment, a well known effect seen 
in simulations of decaying incompressible MHD turbulence 
where the minority species usually decays more rapidly than 
the dominant species causing the ma gnitude of the normalized 
cross-helicity to increase with time 



ignitude or tt 

ei mm. 



In freely decaying turbulence, dynamic alignment occurs 
whenever the total energy decays more rapidly than the cross- 
helicity, that is, e > \s c \, where the cascade rate of cross- 
helicity e c may be positive or negative. From the relations 
e > and e — e ± e c , it follows that dynamic alignment 
occurs if and only if e + > and e~ > 0. If w + > ur , it is 
not necessary that r + > r~ , only that 



r 1 e 
— > 



+ • 



(52) 



as can be seen from equation (f3Tb . Therefore, even though 
< r~ may seem counter-intuitive, it is not 



+ 



the relation r 
inconsistent with dynamic alignment. 



V. SUMMARY AND CONCLUSIONS 

Observations of scale dependent alignment of velocity and 
magnetic field fluctuations 5v± and 5b± in the solar wind 
suggest that this effect must be included in any theory of so- 
lar wind turbulence ITol [Till . Perez and Boldyrev lfl7Tl have 
recently discussed a theory of imbalanced turbulence that 
includes scale dependent alignment of the fluctuations 6v± 
and 5b± in the inertial range. We have extended the Perez - 
Boldyrev theory by including the probabilities p and q which 
solar wind observations indicate are not necessarily equal. 
Operationally, the probabilities p and q may be defined as fol- 
lows. Suppose space is covered by a uniform cartesian grid or 
three dimensional mesh. At each grid-point one may compute 
the fluctuations Sv± and 5b± and the angle between them 9. 
If the angle lies in the range < 6 < tt/2, then the fluctuation 
is positively aligned and if tt/2 < 9 < tt, then the fluctuation 
is negatively aligned. By counting the number of positively 
and negatively aligned fluctuations in a large volume V, much 



larger than the lengthscales of the turbulent eddies, the prob- 
abilities p and q may be defined as the fractional numbers of 
positively and negatively aligned fluctuations in the volume 
V. 

The phenomenological theory developed in this paper was 
guided primarily by two new solar wind observations. It 
should be noted that both of these solar wind observations are 
necessary for the development of the theory. At first glance, it 
may seem that the condition a c = const implies that p and q 
are both constant. Or that these two conditions are somehow 
equivalent. However, the relation (w + /w~) 2 ~ pvi/qv 2 , 
equation (|24| >. shows that p/q can vary with the lengthscale 
even if w + /w~ is constant. Therefore, it is essential to have 
separate observations of the scale invariance of a c and the 
scale invariance of p and q to support the theoretical frame- 
work developed here. 

In summary, using estimates of the cascade times derived 
from the nonlinear terms in the incompressible MHD equa- 
tions and two new observational constraints derived from 
studies of solar wind data, we have constructed a generaliza- 
tion of Boldyrev's theory [7] that depends on the three param- 
eters w + /w~, e + /e~, and p/q. The theory reduces to the 
original theory of Boldyrev O] when w + = w~, e + = e~, 
and p — q since in this limit 9^ = 0j_ and the cascade times 
d32l become equal to those of Boldyrev 0]. For imbalanced 
turbulence w + ^ ur , p ^ q, and the theory predicts the 
scaling laws w ± oc A^ 4 , oc X 1 / 4 , and oc X 1 / 2 - 
Interestingly, the scaling laws for balanced and imbalanced 
turbulence are the same. The perpendicular energy spectrum 
defined by k±E ± ~ |w ± | 2 has the inertial range scaling 
E ± oc fc~ 3/2 with 



El 
E~ 



l-a c 



const. 



(53) 



The theory assumes that the cascades for positively and neg- 
atively aligned fluctuations are both in a state of critical bal- 
ance d42b . although they are governed by different timescales, 
and that the eddy geometry is scale invariant. The positively 
aligned fluctuations occupy a fractional volume p and the neg- 
atively aligned fluctuations occupy a fractional volume q so 
that the energy cascade rate is 



or, equivalently, 



u*. 

)— 



4 T H 



turf 

4t- 



(54) 



(55) 



In the discussion following equation (l35l l it was shown that 
at most three of the ratios w + /w~, e + /e~, and p/q can be in- 
dependent. However, the two ratios w + /w~ and e + /e~ can- 
not be independent since in the case of homogeneous steady- 
state turbulence w + — w~ implies e + = e~ and vice versa. 



8 



This is because the injection of cross-helicity into the system, 
e c 7^ or e + ^ s~, will create a nonzero cross-helicity spec- 
trum and a cascade of cross-helicity from large to small scales 
which implies a net accumulation of cross-helicity within the 
volume (cr c ^ 0). Hence, at most two of the ratios and 
w + /w~ and p/q are independent. Whether p/q can be ex- 
pressed in terms of w + /w~ and e + /e~ is an open question. 
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